Magnetic heat conductivity in CaCu 2 03: linear temperature dependence 
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We present experimental results for the thermal conductivity k of the pseudo 2-leg ladder material 
CaCu2 03. The strong buckling of the ladder rungs renders this material a good approximation to a 
5=1/2 Heisenberg-chain. Despite a strong suppression of the thermal conductivity of this material 
in all crystal directions due to inherent disorder, we find a dominant magnetic contribution K mag 
along the chain direction. K mag is linear in temperature, resembling the low-temperature limit of the 
thermal Drude weight Dth of the S = 1/2 Heisenberg chain. The comparison of K mag and Dth yields 
a magnetic mean free path of Z mag « 22 ± 5 A, in good agreement with magnetic measurements. 

PACS numbers: 75.40.Gb,75.10.Pq,66.70.+f,68.65.-k 



Recently, the magnetic heat transport of low- 
dimensional quantum spin systems has become the focus 



of numerous studies [l|. 
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since intriguing prop- 
erties have been found. Firstly, a substantial magnetic 
contribution to the thermal conductivity k in addition to 
a regular phononic background has been experimentally 
established for spin chain and ladder compounds as well 
as two-dimensional (2D) antiferromagnets such as the in- 
sulating parent compounds of superconducting cuprates. 
It can thus be considered a generic feature of quasi low- 
dimensional magnetic materials. Secondly, the surpris- 
ingly large magnetic contribution to k of spin ladder 
materials observed in the case of (Sr, Ca, La)i4Cu2404i 
0, 0- IE 0] has triggered extensive theoretical work on 
possible ballistic heat transport in spin chains and lad- 
ders [3, EI EE EJ EE EE dl . 



The analysis of experimental data for n can be quite 
involved, in particular when phonon and magnetic en- 
ergy scales do not separate as is the case for the spin 
chain compounds Sr2Cu03 and SrCu02 0,0. The sit- 
uation in the cases of spin ladder materials and the two- 
dimensional cuprates is quite fortunate: here, a clear sep- 
aration of phonon (ft p h) and magnetic (ft mag ) contribu- 
tions allows for a robust determination of K mag itself. It 
is important to note that K ma „ of the spin ladder com- 
pounds (Sr, Ca, La)i4Cu2404i [l|, H, H, 3 shows thermally 
activated behavior at low temperature T dominated by 
the large spin gap of two- leg ladders [l], , while in the 
case of La2Cu04 Q K mag is found to be proportional 
to T 2 , which is the leading contribution in T to the spe- 
cific heat of a 2D square lattice antiferromagnet. Thus in 
these two cases, the experimentally observed K mag clearly 
exhibits intrinsic properties of the underlying spin mod- 
els. 

Theoretically, a consistent picture for thermal trans- 



port in spin-1/2 Heisenberg chains has only recently 
emerged: the integrability of this model results in a diver- 
gent K mag [LJ, |l5|, |l_6( which is described by the so-called 
thermal Drude weight D t h multiplied by a delta function 
at zero frequency. At low T this D t h depends linearly on 
T [IH and such behavior is more generally expected for k 
of any spin chain model with gapless excitations 

nana. 

It is the purpose of this Letter to present the first 
experimental example of a quasi one-dimensional (ID) 
quantum magnet (namely CaC^Oa) which exhibits a 
linear T-depcndcnce of K mag over a wide T-range. This is 
an eye-catching result as it resembles the intrinsic trans- 
port properties of a spin chain. The extraction of K mag 
from the experimental data is very accurate since K p h 
of this material is strongly suppressed due to inherent 
disorder. We utilize well-known expressions for D t h to 
extract information on the scattering processes. Our 
analysis yields a T-independent magnetic mean free path 
( mag w 22 ±5 A. The mean separation of magnetic defects 
along the chain direction as determined from magnetiza- 
tion measurements is of the same order of magnitude. 

The space group of CaCu203 is Pmmn with lattice 
constants a = 9.946 A, b = 4.079 A, and c = 3.460 A 
24j. The structure basically consists of CU2O3 structural 
units arranged in the afc-plane with the geometrical form 
of buckled two-leg ladders running along the 6-direction. 
The Cu2C>3-planes are stacked along the c-direction and 
are separated from each other by layers of Ca-ions. The 
nearest neighbor magnetic exchange coupling of the Cu 2+ 
spins along the 180° Cu-O-Cu bonds which form the lad- 
der legs in the 6-direction is large and has been estimated 
as J/ks ~ 2000 K [25|. The intra-ladder magnetic rung- 
coupling along the buckled Cu-O-Cu bonds (123° bond- 
ing angle) in the a-direction is much smaller and believed 
to be in the range J±/k B « 100-300 K [H, El . The 
inter-ladder coupling along the c-axis has been estimated 
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FIG. 1: K a (□), Kf, (O) and k c (A) of CaCu203 as a function 
of T. The dashed and solid lines represent a linear fit of the 
experimental data in the range 100-300 K and the estimated 
K ma g in this range. Extrapolations of re mag towards low T 
(assuming a T-independent Z mag as extracted for T > 100 K) 
corresponding to a finite (A = 3meV) and a vanishing spin 
gap are represented by dotted and dashed-dotted lines, re- 
spectively. Inset: Kb from two different measurements (open 
symbols depict the same curve as in the main panel). 



to be of the same order of magnitude [HHl- Along the 
a-axis a weak (< 100 K) and frustrated nearest neighbor 
inter-laddcr coupling is expected [25| and an even weaker 
(8 — 30 K) second nearest neighbor inter-laddcr coupling 
has been suggested to be mediated via excess Cu 2+ ions 
on interstitial positions [2^ |. Despite CaCu2C>3 being a 
ladder-like material, where in principle a gapped non- 
magnetic ground state is expected [28|, the material or- 
ders antiferromagnetically at T/v ~ 25 K [25, 26|, 2^] and 
inelastic neutron scattering reveals a spin chain-like ex- 
citation spectrum with the upper bound for a spin gap 
A w 3 mcV (32}- We therefore consider CaCu203 as an 
ensemble of weakly coupled spin chains. 

We have grown single crystalline CaCu2C>3 by the 
traveling solvent floating zone method [30l |. Samples 
with typical dimensions of around 3.5 mm length along 
the thermal current direction and 1 mm 2 cross section 
were cut from a crystal and the thermal conductivity 
of CaCu203 has been measured along the a, b and in- 
direction (re a , ret, and re c ) as a function of T in the range 
7-300 K. We used a standard four probe technique where 
errors due to radiation loss are minimized. Magnetization 
has been measured using a superconducting quantum in- 
terference device (Quantum Design MPMS XL5). 

Fig-fflshows K a , Kb and k c of CaCu2C>3 as a function of 
T. Since the material is insulating, electronic heat con- 
duction is negligible and we therefore expect these com- 
ponents to originate from phononic heat conduction plus 



a possible magnetic contribution. The thermal conduc- 
tivities perpendicular to the chain direction (re and k c ) 
only exhibit a weak T-dcpcndencc and share absolute val- 
ues K 4 Wm _1 K _1 ) which arc 1-2 orders of magnitude 
smaller than re p h of other chemically undoped cuprates 
such as SrCuC>2 Q or Sri 4 Cu2404i Q. Instead of a pro- 
nounced phononic low-T peak, which is usually found in 
such cases, only a small peak is present in k c while no 
peak is found in K a . In fact, such strongly suppressed 
k is typical for re p h with a high phonon scattering rate 
[3H . Indeed, substantial phonon-defect scattering must 
be present in CaCu2C>3 due to inherent structural disor- 
der induced by a significant Ca and oxygen deficiency be- 
ing balanced by excess Cu 24 , 2?J ■ We therefore consider 
k perpendicular to the chains to be purely phononic. 

A completely different behavior is observed for re par- 
allel to the chains. At low T < 50 K Kb resembles re c 
although it is about three times larger. However, Kb 
steeply increases at T > 50 K with this increase being 
linear in T for T > 100 K. Such a strong increase of re 
with rising T cannot be understood in terms of conven- 
tional phonon heat conduction by acoustic phonons. Dis- 
persive optical phonons could in prin ciple give rise to an 
increase of re p h with increasing T [32j and possibly play a 
role in the T dependence of re a . Heat transport by opti- 
cal phonons is however unable to account for the rather 
large re& at 300 K since this would require unrealistically 
large phonon mean free paths [38j |. In analogy with ob- 
servations in other cuprate materials 0, 0, d, 0, @, 0, @] , 
we therefore conclude that the strong increase of Kb and 
the resulting anisotropy of the re tensor originate from ID 
heat transport due to magnetic excitations in the weakly 
coupled spin chains. While these observations and the 
following analysis represent the central results of this pa- 
per, we wish to point out an additional feature of the 
thermal conductivity. As depicted in the inset of Fig. 1, 
we observe variations of re as measured during different 
runs which are beyond the statistical error of the typi- 
cal measurement. Similar has been observed in other ID 
quantum magnets as well [ll , 3^, 34 , 35[ and remains an 
open issue which merits further investigation. 

In order to separate the phononic (re p h,&) and magnetic 
(^mag) parts of Kb we assume that re p h,& is of a similar 
magnitude and exhibits a similar T-dependence to the 
purely phononic re a and re c . Possible errors due to the 
crudeness of this estimation become small at T > 100 K 
where the strongly increasing magnetic part of Kb be- 
comes clearly larger than any possible phononic ther- 
mal conductivity. Since re a and re c are only weakly T- 
dependent and Kb increases linearly in this T-regime, it 
is natural to conjecture re p h,b ps const and re mag oc T 
at T > 100 K. Indeed, a linear fit in the range 100- 
300 K (dashed line in Fig. Q]) describes the data al- 
most perfectly and yields re p h,& = 1.2 Wm _1 K _1 and 
K raag = 0.055 Wm-'K" 2 x T. We plot the extracted 



for T > 100 K as a solid line 
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at lower T cannot be inferred from our data. Return- 
ing to the inset we note that the alternative curve for 
Kf, also increases linearly with T for T > 100 K yielding 
K ph ,6 = 2.2 Wm^K" 1 and K mag = 0.041 Wm^K^xT. 
We therefore assume that the difference in the two results 
for Kb represents the error in determining the intrinsic k 
of our sample. fc p h,& could certainly exhibit a slight In- 
dependence, either increasing (as K a ) or decreasing (as 
normally expected for K p h) with rising T. These uncer- 
tainties are irrelevant to our further analysis. 

Qualitatively, the linear increase of K mag with T is a re- 
markable result as it directly reflects the low-temperature 
behavior of the thermal Drude weight of a Heiscnbcrg 
chain 15, 1^1, where (in SI units) 
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As the exchange integral J/Hb along the chain is of the 
order of 2000 K, we may safely neglect any deviations 
from the linear behavior that become relevant at tem- 
peratures T > 0. 15 J/k B [HEil, i- e -> T ~ 300 K - 

As mentioned earlier, the structure of CaCu203 sug- 
gests a model of weakly coupled chains. We have checked 
that any deviations of R mag from a linear T-dependence 
due to the presence of a small spin gap can only occur 
in the low temperature regime which is dominated by 
phonons and thus experimentally hardly observable. Us- 
ing a Boltzmann-type expression for K mag (see Eq.(l) in 
Ref. 0) to estimate the effect of a small gap, we fur- 
ther find that at T > 100 K a pure spin chain and a 
weakly coupled ladder for J±/J < 0.05 both result in 
the same linear T-dependence (cf. Fig. []}. For the lat- 
ter, A < 3 meV since A w 0.4J_i_ [36$ • It is thus justified 
to use Eq. ((TJ) to analyze our data at high T. 

In the presence of external scattering, the thermal con- 
ductivity of a single chain /t mag is rendered finite with 
a width ~ 1/t and may be approximated by fc mag = 
Dth t/tt. Combining this with Eq. JT]) to calculate the 
magnetic mean free path ' mag = vt from our experimen- 
tal data for re r 



yields 



3 ft 



'mag 



7T k 2 B N 



T 



(2) 



where N = A/ac is the number of spin chains per unit 
area in the crystal. A kinetic approach as described in 
Ref. yields the same result in the low temperature limit. 
Note further that in our model, r is an energy indepen- 
dent quantity which seems to be a reasonable approxi- 
mation, since even at T = 300 K <C J/ks only spinons 
in a small region of the Brillouin zone can contribute to 
K mag . In this region of the Brillouin zone v w Jair /2ft is 
constant and hence ' mag can also be regarded as indepen- 
dent of energy. From Eq. [2] and the experimental data 
we obtain i mag = 22 ± 5 A for the entire range 100-300 K 
[39l ] , corresponding to about 5 — 6 lattice spacings. 
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FIG. 2: Magnetization of CaCu203 as a function of tempera- 
ture with a magnetic field fioH = 1 Tesla parallel to the three 
crystallographic axes. The solid lines represent Curie- Weiss- 
type fits to the data in the range 100-360 K which yield ~ 3% 
free spins with respect to Cu. Inset: Inverse magnetization 
after subtracting a constant Mo which roughly accounts for 
van Vleck and chain magnetism [25l |. 



Within the framework of Boltzmann-type kinetic mod- 
els, T-independent mean free paths for magnetic excita- 
tions have already been observed in the spin ladder mate- 
rial (Sr, Ca,La)i 4 Cu 2 40 4 i 0,0, HQ and the 2D antifer- 
romagnet La2Cu04 Q where in the low-T regime scat- 
tering off defects dominates over other possible scattering 
mechanisms. In analogy with these cases it seems rea- 
sonable to explain our result of a T-independent Z mag in 
CaCu203 by dominant spinon-defect scattering as well. 
However, a constant ' mag over such a large T-range is 
very surprising because spinon-phonon scattering should 
become increasingly important at higher T and eventu- 
ally lead to a T-dependent mean free path. Despite this 
intuition the data suggest that at 300 K the probability of 
a spinon scattering off a phonon is still much lower than 
that of scattering off a defect. It is therefore natural 
to conclude that the density of relevant defects is very 
high in the material, resulting in scattering off defects 
being much more important than other scattering pro- 
cesses. This is consistent with our quantitative result for 
^mag which is about 1-2 orders of magnitude smaller than 
in spin chain systems with similar magnetic exchange, 
e.g. SrCu02 and Si^CuOa, where much larger and T- 
dependent mean free paths have been found [y, 0] • 

It is obvious to search for the origin of the high den- 
sity of static scattering centers as suggested by our data 
in the strong off-stoichiomctry of CaCu203. The oxygen 
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deficiency is clearly most relevant for creating vacancies 
within the CU2O3 chain structures which must create lo- 
cal scattering sites. In order to estimate the density of 
such scattering sites we consider the T-dependence of the 
magnetization M of our sample, shown in Fig. [5J which 
is in good agreement with previous results for crystals 
from a different sample growth 2^, 26[. M(T) exhibits a 
sharp anomaly at the Neel temperature TV ~ 25 K and 
a Curie-like behavior at T > T/v indicating ~ 3% of free 
spins (with respect to Cu) in the material. Kiryukhin et 
al. suggested |25j that these free moments are located 
directly in the chains and originate from chain interrup- 
tions. According to their analysis the high-T Curie tail 
in the magnetization originates from chain segments with 
a length between ~ 40 and 80 A. In contrast, a recent 
electron spin resonance study by Goiran et al. [26| sug- 
gests that the free moments are more likely to arise from 
excess Cu 2+ ions where each of these ions resides on an 
interstitial site in the vicinity of an oxygen vacancy in a 
neighboring chain structure. Within the latter scenario 
a lower limit for the density (per unit length along the 
^-direction) of oxygen vacancies n > 0.03/6 can be in- 
ferred from the M(T) data. This yields an upper limit 
for the mean distance d between the induced local scat- 
tering sites within a chain of d = ^ < 68 A. In this 
estimate we take into account that each vacancy creates 
a local structural distortion which affects at least two 
chains since the individual chain structures are strongly 
interwoven. In either case the extracted mean length of 
non-distorted chain segments d should be regarded as an 
upper limit for the magnetic mean free path since not ev- 
ery chain distortion affecting K mag necessarily contributes 
to the magnetization. This agrees reasonably well with 
the analysis of the thermal conductivity as / mag is found 
to be somewhat smaller (by a factor 2-3) than d. 

We mention that the T-independent scattering rate as 
found in our experiment is in conflict with recent calcula- 
tions of a scattering rate oc T^ 1 for impurities which 
induce slight disorder in the magnetic exchange coupling 
j22l |. A possible reason for this discrepancy could be re- 
lated to the nature of the defects in our sample. Very 
likely the degree of the distortion of the chain at the 
actual scattering sites is too large for this model to be 
applicable. Detailed experimental and theoretical inves- 
tigations are underway to clarify this issue. 

To conclude, we have studied the thermal conductivity 
of CaCu203 as a function of temperature T. The thermal 
conductivity parallel to the chains exhibits a pronounced 
linear increase with increasing T which we attribute to 
magnetic heat transport within the chains of the mate- 
rial. The linear increase resembles the intrinsic thermal 
conductivity of a spin chain with a constant scattering 
rate. We extract a value for the magnetic mean free path 
i mag « 22 ± 5 A which is in reasonably good agreement 
with the mean distance between magnetic defects in the 
material as determined from magnetization data. 
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